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We discuss a classical nonlinear oscillator, which is proved to be a superin- 
tegrable system for which the bounded motions are quasiperiodic oscillations 
and the unbounded (scattering) motions are represented by hyperbolic func- 
tions. This oscillator can be seen as a position-dependent mass system and we 
show a natural quantization prescription admitting a factorization with shape 
invariance for the n = 1 case, and then the energy spectrum is found. Other 
isochronous systems which can also be considered as a generalization of the 
harmonic oscillator and admit a nonstandard Lagrangian description are also 
discussed. 
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The harmonic oscillator is a system playing a privileged role both in classical 
and quantum mechanics. It is almost ubiquitous in Physics and appears in many 
physical applications running from condensed matter to semiconductors (see e.g. 
for references to such problems). The dynamical evolution of the classical 
^ . system in one dimension is given by 

P. , dq _ dv _ 2 

dt ' dt ^ ' 

and admits a Lagrangian formulation with L = (1/2) [v 2 — o; 2 g 2 ), the general 
solution of the equations of motion being 

Q = Qo cos ujt — — sinujt = Acos(ujt + tp) 

UJ 

and therefore the solutions are periodic with angular frequency u, while A and cp 
are arbitrary. This is the main characteristic of the classical system. As far as the 
quantum system is concerned, the eigenvalues of the Hamiltonian, which is given 
by H = (1/2) (p 2 +oj 2 q 2 ), are equally spaced. We should also remark that the 
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natural extensions to two dimensions, given by 

H = \{ P l+V 2 y ) + \(^x 2 + uW) 

admits two constants of motion in involution, 1\ = E x = \ (p 2 + u> 2 x 2 ), I2 = E y = 
I {p 2 + uj\ y 2 ), and therefore it is completely integrable in the sense of Liouville. 
Moreover it has been proved that, when u\ and 102 are rationally related, i.e. 
uj\ = niwo, W2 = ri2 Uq, with m, n-z G N, there exist a new constant of motion and 
the system is superintegrable. Actually the complex function, J = K™ 2 (K*) ni 
with K x = p x + iniujQX and K y = p y + % ri2 ujq y, is a constant of the motion. 

Our aim is to comment on some possible generalizations of this system from the 
perspective of the theory of the symmetry, i.e. trying to preserve the fundamental 
symmetry properties. 



2 A position-dependent mass nonlinear oscillator 

A often used generalization was proposed by Mathews and Lakshmanan [2J |3| as 
a one-dimensional analogue of some models of quantum field theory [IJ [S] . It is 
described by a Lagrangian 

1 / 1 ^ - -„2 „2 2\ 



which can be considered as an oscillator with a position-dependent effective mass 
m = (1 + Ax 2 )^ 1 (see e.g. [SUZ] and references therein). It was proved that the 
general solution is also q{t) = Acos{ujt + ip), but now the amplitude A depends on 
the frequency. More explicitly uj 2 (1 + A A 2 ) = a 2 . Note also that this Lagrangian 
is of mechanical type, the kinetic term being invariant under the tangent lift of 
the vector field 

X X (X) = Vl + Ax 2 |- . 

It was recently shown in jS] that there is a generalization to n dimensions 
preserving the symmetry characteristics. In particular the two-dimensional gen- 
eralization studied in [H] was given by the Lagrangian 

1 / 1 



L{\) = 2 (r+A7V [< + v y + \(xv y -yv x ) 2 -a 2 r 2 \ , r 2 = x 2 + y 2 , (2) 

and it was shown to be not only integrable but also superintegrable. This is 
the only generalization to n dimensions for which the kinetic term is a quadratic 
function in the velocities that is invariant under rotations and under the two 
vector fields generalising the symmetries of the one-dimensional model, i.e. 

d / — d 



Ax(A) = Vl + Ar 2 — , A 2 (A) = ^1 + Ar 2 — . 
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This is valid for any value of A. However, when A < 0, A = — |A|, this function 
has a singularity at 1 — |A| r 2 = and we restrict our dynamics to the interior of 
the circle x 2 + y 2 < 1/|A|. 

These two vector fields close with the generator of rotations, 

d d 

X J ' = X TT ~y-fT ' 
ay ox 

on a Lie algebra 

[X 1 (X),X 2 (X)] = XXj, [X 1 (X),X J ] = X 2 (X), [X 2 (X),X J ] = -X 1 (X). 

which is isomorphic either to ^O^R), when A > 0, or to SO{2, 1), when A < 0, 
or finally to the Euclidean group in two dimensions when A = 0. 

The important property shown in [5j is that this bidimensional nonlinear har- 
monic oscillator is completely integrable, because one can show that, if K\ and 
K 2 are the functions 

K 1 = P 1 {X)+ia - X = , K 2 = P 2 (X) + ia 



Vl + Xr 2 ' v ' Vl + Ar 2 ' 

with 

r>t\\ v x -XJy Vy + XJx 

Pl ^ = ^r^T> P2 ^ = n=^n> J = xv y -yv x , 
V1+ Ar v 1+Ar 

then the complex functions Kij defined as K^j = KiKj, i,j = 1,2, are constants 
of motion. In fact the time-evolution of the functions K\ and K 2 is 

d K ia K d K - ia K 

from which we see that the complex functions are constants of the motion. 
Therefore the system is superintegrable with the following first integrals of motion 

Ji(A) = | K x | 2 , I 2 {X) = \K 2 \ 2 , h = %K l2 ) = a (xv y - yv x ) . 

The Legendre transformation for a two-dimensional Lagrangian system of me- 
chanical type with kinetic term as in © is given by 

(1 + A y 2 )v x - X xyvy _ (1 + A x 2 )v y - X xyv x 



1 + Ar 2 ' iy 1 + Ar 2 

(note that xp y — yp x = xv y — yv x ) and the general expression for the corresponding 
A-dependent Hamiltonian is 

H(X) = \ [pl+p 2 y + ^(xp x + yp y ) 2 ] +^a 2 V(x,y), (3) 
and hence the associated Hamilton-Jacobi equation is 

fds\ 2 fds\ 2 x / ds ds\ 2 2Tr , , 
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This equation is not separable in (x, y) coordinates, but it was shown in [H] 
that there exist three particular systems of orthogonal coordinates, and three 
particular families of associated potentials, in which such Hamiltonians admit a 
Hamilton-Jacobi separability. The first system of coordinates is given by 



\%x > y) i %x 



(5) 



for which the Hamilton-Jacobi equation becomes: 

(1 + \zl) (g) 2 + (1 + Xy 2 ) 2 (H) 2 + a 2 (1 + Xy 2 )V = 2(1 + Xy 2 )E , 

and therefore the Hamilton-Jacobi equation is separable if the potential V(x, y) 
can be written on the form 



V 



l + Xy 2 



+ W 2 (y) 



(6) 



The potential is therefore integrable with the following two quadratic integrals of 
motion 



I 2 (X) = (l + Xr 2 )p 2 -XJ 2 + a 2 (w 2 (y) 



Xy 2 



l + Xy 2 



In a similar way, one can see, using coordinates (x, z y ) with z y = y(l+Ax 2 ) -1 / 2 , 
that the Hamilton-Jacobi equation is separable when the potential V(x,y) is of 
the form 



V = Wi(x) + 



l + Xx 2 ' 

and the potential is integrable with the following two quadratic first integrals: 

Ax 2 



(7) 



h(X) 
h(X) 



(1 + A r 2 )p 2 x - X J' 2 + a 2 ( Wi (x 



1 + Ax< 



W^Zy) 



(l + Xr 2 )p 2 + a 2 W 2 {zy). 
Finally using polar coordinates (r, eft) the Hamiltonian H(\) is written 



H(X) 



(1 + Xr 2 )p 2 + P 4 



(8) 



and the Hamilton-Jacobi equation is given by 



(1 + Ar 2 



(dS\ 2 1 fdS^ 2 



- 



\ dr 



+ 



+ a 2 V(r,4)) =2E. 



Two important examples of nonlinear oscillator 



5 



Then the equation admits separability when the potential V is of the form 

V = F(r) + G((/))/r 2 . (9) 
Such a potential V is integrable with the following two quadratic first integrals: 

ii(A) = (l + Xr 2 )p 2 r + 

J 2 (A) = 4 + a 2 G(0). 
Consequently, the potential 



1 - r 2 



1-r 2 



x 2 + y 2 



2 Vl + A(x 2 +y 2 

is super-separable since it is separable in three different systems of coordinates 
(z x ,y), (x,z y ), and (r,4>) because 



V- 



a 



1 



2 Vl + Ay 2 / Ll + Az 2 



zl 



a 



1 



2 Vl + Ax' 



x 2 + 



1 + Az 2 



a 



2 Vl + Ar 2 



3 The one-dimensional quantum nonlinear oscillator 

We now consider the quantum case and restrict ourselves to the one-dimensional 
case. The first problem is to define the quantum operator describing the Hamilto- 
nian of this position-dependent mass system, because the mass function and the 
momentum P do not commute and this fact gives rise to an ambiguity in the or- 
dering of factors. It has recently been proposed to avoid the problem by modifying 
the Hilbert space of functions describing the system [2]. More explicitly we can 
consider the measure dfj, = (1 + Ax 2 ) _1//2 dx, which is invariant under the vector 
field X X (X) = V 1 + A x 2 d/dx, for then the operator P = — i V 1 + A x 2 d/dx is 
selfajdoint in the space L 2 (R,dfi). In the case of the nonlinear oscillator in which 
we are interested we can consider the Hamiltonian operator 



1 . , 2 N d 2 1 d 1 a 2 x 2 
-2 (1 + AX) ^-2 AX ^ + 21TA^- 



(10) 



The spectral problem of such operator can be solved by means of algebraic tech- 
niques. We first remark that if f3 is such that a 2 = f3((3 + A), then H[ = H\ — (3/2 
can be factorized product H[ = A\(3) A(f3) and 



a 1 ( r. — ; — o d f3x 
A = —F Vl + Ax 2 — + . H = 



for which its adjoint operator is 



At 



1 



V2 V 



f-vT 



+ Ax 2 



d 



+ 



fix 



dx Vl + Ax 2 



(11) 



(12) 
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The important point is that the partner Hamiltonian H' 2 = A(j3) A* (/3) is related 
to H[ by H' 2 (P) = HUPi) + R(fix) with f3 1 = f((3) and where / and R are the 
functions /(/?) = (3 — A and R(/3) = (3 + (1/2). Hamiltonians admitting such 
factorization ^Oj and related to its s partner in such a way are said to be shape 
invariant and their spectra and the corresponding eigenvectors can be found by 
using the method proposed by Gendenshtein ^2 E] (see also ^H] for a modern 
presentation based on the Riccati equation). Therefore, as the quantum nonlinear 
oscillator is shape invariant, we can develop the method proposed in (111 I12| for 
finding both the spectrum and the corresponding eigenvectors. The spectrum is 
given by 

n 2 1 
E n = n(3-—\ + -(3. 

The existence of a finite or infinite number of bound states depends up on the 
sign of A as also discussed in [S] . 



4 Periodic motions and another nonlinear oscillator 

Another possible generalization of the harmonic oscillator would be to look for 
alternative isochronous systems. For instance one can consider a potential 



U{x) 



Ui{x) ifx<0 
U 2 (x) if x > , 



where U^ix) is an increasing function and U\(x) is a decreasing function, and try 
to determine the explicit functions U\ and U2 in order to have an isochronous 
system. The problem of the determination of the potential when the period is 
known as a function of the energy was solved by Abel Jl] . When the potential is 
symmetric the solution is unique. Therefore the only symmetric potential giving 
rise to isochronous motions around the origin is the harmonic oscillator. The 
isotonic oscillator is also symmetric and isochronous, but the origin is a singular 
point and not a minimum of the potential. Other nonsymmetric potentials can 
be used, for instance a potential given by 

U\(x) = uj\ x 2 , U2 (x) = to 2 x 2 . 

If we want to find more general solutions for the symmetric case we may con- 
sider Lagrangians of a nonstandard mechanical type, in which there is no potential 
term. These more general Lagrangians can also be relevant in other problems. 
For instance another interesting oscillator-like system has recently been studied 
by Chandrasekar et al ^Sj. As mentioned in that paper the oscillator-like system 
admits a Lagrangian formulation. We recall that there are systems admitting a 
Lagrangian formulation of a nonmechanical type. As an example we can consider 
Q = R as the configuration space and the Lagrangian function JH] 

L(x,v) = (a(x)v + U(x)y 1 , (13) 
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which is singular in the zero level set of the function (p(x, v) = a{x) v + U{x). 
Then the Euler-Lagrange equation is 

2a(x) [a' (x)v 2 + U'(x)v + a{x)a] 
a (x)v + U (x) — a (x) v = 



a(x)v + U(x) 

where v and a denote the velocity and the acceleration, respectively. 

This is a conservative system the equation of motion of which can be rewritten 

as 

3 1 

[a(x)} 2 x + a(x) a'(x) x 2 + - a(x) U'(x) x + - U(x) U'(x) = . 

The energy is given by Ei(x,v) = —[2a(x)v + U(x)][a(x) v + U(x)]~ 2 . In 
particular, when a(x) = 1, the Lagrangian is L(x,v) = [v + U(x)]~ 1 and the 
Euler-Lagrange equation reduces to 

x + -U'(x)x+-U(x) U'(x) = (14) 

and the energy function turns out to be El(x,v) = — [2v + U(x)][v + U(x)]~ 2 . 
When U(x) = kx 2 , the equation is 

x + 3kxx + k 2 x 3 = , 

and the energy is El = — [2v + fcx 2 ][« + kx 2 ]^ 1 . It can be seen from the energy 
conservation law that the general solution is 

2t 



kt 2 - E 

The two-dimensional system described by L(x,y,v x ,v y ) = [v x + kix 2 ]^ 1 + 
[v y + k2 y 2 ] -1 is sup erintegr able. Actually not only the energies of each degree of 
freedom are conserved but also the functions |16j 

T x y T &2 ki k\k2xy 



v x + kix 2 v y + k 2 y 2 ' v x + kix 2 v y + k 2 y 2 {v x + h x 2 )(v y + k 2 y 2 ) 

Another example is that of a nonlinear oscillator for which we were looking. 
The following Lagrangian depending on the parameter u 

L(x,v;uj) = - I — f , (15) 

kv x -\- k x -\- w z 

produces the nonlinear Euler-Lagrange equation 

x + 3kxx + k 2 x 2 + uj 2 x = , 

which is the nonlinear oscillator system recently studied by Chandrasekar et al 
[15] . and the energy is Er = -[2 k v x +k 2 x 2 +oj 2 ) [(k v x +k 2 x 2 +uj 2 ))~ 2 . The general 
solution for the dynamics, which can be found from the energy conservation, is 



uj\ E sin(u;i + (f>) 
1 — k \[E cos(u;£ + 4>) 
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We have recently been able to prove |16j that in the rational case of the two- 
dimensional problem, for which lo\ = n\ loq and 0J2 = ^2^0, the system is superin- 
tegrable as it was the case for the harmonic oscillator. To introduce the additional 
constants of motion we define 

v x + k\ x 2 + i ni loq x v y + k 2 y 2 + in 2 uJ y ,, c s 

= 7n ^ 9 9 , llV2 = "J To o 9 9 5 

k\ v x + kfx z + n\ Wq k 2 v y + k^y + ^2 w o 

and then the complex function K™ 2 (K^)™ 1 is a constant of the motion. 

In summary, not only position-dependent mass generalizations of the harmonic 
oscillator can be interesting but there exist also systems described by Lagrangians 
of non-mechanical type which preserve the property of superintegrability for the 
harmonic oscillator with rationally related frequencies. This example points out 
the importance of the study of such non-standard Lagrangians. 
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